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Abstract. This article summarizes theoretical predictions for the density and isospin dependence of the 
nuclear mean field and the corresponding nuclear equation of state. We compare predictions from micro- 
scopic and phenomenological approaches. An application to heavy ion reactions requires to incorporate 
these forces into the framework of dynamical transport models. Constraints on the nuclear equation of 
state derived from finite nuclei and from heavy ion reactions are discussed. 
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1 Introduction 

Heavy ion reactions provide the only possibility to reach 
nuclear matter densities beyond saturation density po — 

0. 16 fm -3 . Transport calculations indicate that in the low 
and intermediate energy range E\ a h ~ 0.1 -r 1 AGeV nu- 
clear densities between 2 3po are accessible while the 
highest baryon densities (~ 8po) will probably be reached 
in the energy range of the future GSI facility FAIR be- 
tween 20 -j- 30 AGeV. At even higher incident energies 
transparency sets in and the matter becomes less baryon 
rich due to the dominance of meson production. The isospin 
dependence of the nuclear forces which is at present only 
little constrained by data will be explored by the forth- 
coming radioactive beam facilities at FAIR/GSI 1 , SPI- 
RAL2/GANIL and RIA 0. Since the knowledge of the 
nuclear equation-of-state (EOS) at supra-normal densities 
and extreme isospin is essential for our understanding of 
the nuclear forces as well as for astrophysical purposes, 
the determination of the EOS was already one of the pri- 
mary goals when first relativistic heavy ion beams started 
to operate in the beginning of the 80ties In the follow- 
ing we will briefly discuss the knowledge on the nuclear 
EOS from a theoretical point of view, then turn to the re- 
alization within transport models, and finally give a short 
review on possible observables from heavy ion reactions 
to constrain the EOS. 

2 Models for the nuclear EOS 

Models which make predictions on the nuclear EOS can 
roughly be divided into three classes: 

1. Phenomenological density functionals: These are 
models based on effective density dependent interac- 
tions such as Gogny |4I5| or Skyrme forces (617) or rel- 
ativistic mean field (RMF) models jHj- The number of 



parameters which are fine tuned to the nuclear chart 
is usually larger than six and less than 15. This type 
of models allows the most precise description of finite 
nuclei properties. 

2. Effective field theory approaches: Models where 
the effective interactions is determined within the spirit 
of effective field theory (EFT) became recently more 
and more popular. Such approaches lead to a more 
systematic expansion of the EOS in powers of den- 
sity, respectively the Fermi momentum fcp. They can 
be based on density functional theory (9110) or e.g. on 
chiral perturbation theory [11112113] . The advantage 
of EFT is the small number of free parameters and 
a correspondingly higher predictive power. However, 
when high precision fits to finite nuclei are intended 
this is presently only possible by the price of fine tuning 
through additional parameters. Then EFT functionals 
are based on approximately the same number of model 
parameters as phenomenological density functionals. 

3. Ab initio approaches: Based on high precision free 
space nucleon-nucleon interactions, the nuclear many- 
body problem is treated microscopically. Predictions 
for the nuclear EOS are parameter free. Examples are 
variational calculations (14115) . Brueckner-Hartree-Fock 
(BHF) (Ifill7ll8ll9| or relativistic Dirac-Brueckner-Hartree- 
Fock (DBHF) I20l2ll22l23l24l25l26| calcu lations an d 
Greens functions Monte-Carlo approaches (27128129) . 



Phenomenological models as well as EFT contain parame- 
ters which have to be fixed by nuclear properties around or 
below saturation density which makes the extrapolation to 
supra-normal densities somewhat questionable. However, 
in the EFT case such an extrapolation is safer due to a sys- 
tematic density expansion. One has nevertheless, to keep 
in mind that EFT approaches are based on low density 
expansions. Many-body calculations, on the other hand, 
have to rely on the summation of relevant diagram classes 
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and are still too involved for systematic applications to 
finite nuclei. 



2.1 Mean field theory 

Among non-relativistic density functionals are Skyrme func- 
tionals most frequently used. The Skyrme interaction con- 
tains an attractive local two-body part and a repulsive 
density dependent two-body interaction which can be mo- 
tivated by local three-body forces. We will not consider 
surface terms which involve gradients as well as spin-orbit 
contributions since they vanish in infinite nuclear matter. 
For a detailed discussion of Skyrme functionals and their 
relation to relativistic mean field (RMF) theory see e.g. 
[7|. The EOS of symmetric nuclear matter, i.e. the bind- 
ing energy per particle has the simple form 
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where the first term in represents the kinetic energy 
of a non-relativistic Fermi gas and the remaining part the 
potential energy. To examine the structure of relativistic 
mean field models it is instructive to consider the sim- 
plest version of a relativistic model, i.e. the olq model 
of Quantum Hadron Dynamics (QHD-I) (HOj- In QHD- 
I the nucleon-nucleon interaction is mediated by the ex- 
change of two effective boson fields which are attributed 
to a scalar a and a vector u) meson. The energy density 
in infinite cold and isospin saturated nuclear matter is in 
mean field approximation given by 

e= ^E FQ +^m* D Q S + ^{r V Q 2 + r S Q 2 s } , (2) 



where the Fermi energy is given by Ep — y/k F + m 
We will denote m* D explicitely as Dirac mass in the fol- 
lowing in order to distinguish it from its non-relativistic 
counterpart. The effective mass absorbs the scalar part of 



the mean field 
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the first two terms in (J2J provide the energy (kinetic plus 
rest mass) of a non-interacting relativistic Fermi gas. 

A genuine feature of all relativistic models is the fact 
that one has to distinguish between the vector density g = 
2k F /3ir 2 and a scalar density gs- The vector density is the 
time-like component of a 4-vector current which spa- 
tial components vanish in the nuclear matter rest frame, 
while gs is a Lorentz scalar. The scalar density shows a 
saturation behavior with increasing vector density which 
is essential for the relativistic saturation mechanism. This 
becomes clear when binding energy E / A = e/ g — Mis 
expanded in powers of the Fermi momentum kp 
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The first term in @ contains the kinetic energy of a non- 
relativistic Fermi gas followed by relativistic corrections 



and the remaining terms are the contributions from the 
mean field. In QHD-I the scalar and vector field strengths 
are given by the coupling constants for the correspond- 
ing mesons I5 = g^/m 2 , and Jy = gt/fnt divided by 
the meson masses. The two parameters Js.v are now fit- 
ted to the saturation point of nuclear matter E/A ~ 
— 16 MeV, go ~ 0.16 fnV 3 which follows from the vol- 
ume part of the Weizacker mass formula. The satura- 
tion mechanism requires that both coupling constants are 
large. This leads automatically to the cancellation of two 
large fields, namely an attractive scalar field £$ = —Tsgs 
and a repulsive vector field Sy = FvQ- As is a typical 
feature of relativistic dynamics the single particle poten- 
tial U = m* D /E*S s ~ Zv (E* = ^Jk 2 + m* 2 ), which is 
of the order of -50 MeV, results from the cancellation of 
scalar and vector fields, each of the order of several hun- 
dred MeV. 

However, with only two parameters QHD-I provides a 
relatively poor description of the saturation point with a 
too large saturation density and a very stiff EOS (K=540 
MeV) . To improve on this higher order corrections in den- 
sity have to be taken into account. This can be done in 
several ways: In the spirit of the original Walecka model 
non-linear meson self-interaction terms have been intro- 
duced into the QHD Lagrangian |8I31| . An alternative are 
relativistic point coupling models where the explicit meson 
exchange picture is abandoned. A Lagrangian of nucleon 
and boson fields with point couplings can be constructed 
in the spirit of EFT and expanded in powers of density 
|9I10| . Finite range effects from meson propagators are re- 
placed by density gradients |9I10| . A third possibility is 
density dependent hadron field theory DDRH 32 33 . In 
DDRH the scalar and vector coupling constants are re- 
placed by density dependent vertex functions rsy(kp). 
The density dependence of these renormalized vertices can 
either be taken from Brueckner calculations thus parame- 
terizing many-body correlations |32I33| or be determined 
phenomcnologically 34 35 . In all cases additional param- 
eters are introduced which allow a description of finite 
nuclei with a precision comparable to the best fits from 
Skyrme functionals. Phenomenological density function- 
als provide high quality fits to the known areas of the nu- 
clear chart. Binding energies and rms-radii are reproduced 
with an average relative error of about ~ 1— 5 %. However, 
when the various models are extrapolated to the unknown 
regions of extreme isospin or to super-heavies predictions 
start to deviate substantially. This demonstrates the lim- 
ited predictive power of these functionals. 



2.2 Effective field theory 

When concepts of effective field theory are applied to nu- 
clear physics problem one has to rely on a separation 
of scales. EFT is based on a perturbative expansion of 
the nucleon-nucleon (NN) interaction or the nuclear mean 
field within power counting schemes. The short-range part 
of the NN interaction requires a non-perturbative treat- 
ment, e.g. within the Brueckner ladder summation. The 
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philosophy behind EFT is to separate short-range corre- 
lations from the long and intermediate range part of the 
NN-interaction. This assumption is motivated by the fact 
that the scale of the short-range correlations, i.e. the hard 
core, is set by the p and w vector mesons masses which 
lie well above the Fermi momentum and the pion mass 
which sets the scale of the long range forces. The density 
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Fig. 1. EOS for symmetric nuclear matter obtained from chiral 
one- and two-pion exchange (case 1, solid curve), by adding 
background fields from QCD sum rules, and finally after fine 
tuning to finite nuclei properties (case 3, dashed curve). Figure 
is taken from |12|. 



functional theory (DFT) formulation of the relativistic nu- 
clear many-body problem is thereby analogous to 
the Kohn Sham approach in DFT. An energy functional 
of scalar and vector densities is constructed which by mini- 
mization gives rise to variational equations that determine 
the ground-state densities. Doing so, one tries to approx- 
imate the exact functional using an expansion in classical 
meson fields and their derivatives, based on the obser- 
vation that the ratios of these quantities to the nuclcon 
mass are small, at least up to moderate density. The ex- 
act energy functional which one tries to derive explicitcly 
when using many-body techniques such as Brueckner or 
variational approaches contains exchange-correlation and 
all other many-body and relativistic effects. The DFT 
interpretation implies that the model parameters fitted 
to nuclei implicitly contain effects of both short-distance 
physics and many-body corrections. 

Recently also concepts of chiral perturbation theory 
(ChPT) have been applied to the nuclear many-body prob- 
lem |12I13| . Doing so, the long and intermediate-range in- 
teractions arc treated cxplicitely within chiral pion-nucleon 
dynamics. This allows an expansion of the energy density 
functional in powers of m n /M or in fcp/M. Like in DFT, 
short-range correlation are not resolved explicitely but 
handled by counter-terms (dimensional regularization) (11 j 
or through a cut-off regularization (T3J- Fig. 2] shows the 
corresponding EOS obtained from chiral one- and two- 
pion exchange between nucleons. In order to account for 
the most striking feature of relativistic dynamics, expressed 
by the existence of the large scalar and vector fields, in 



Refs. |12IR-?| iso-scalar condensate background nucleon self- 
energies derived from QCD sum rules have been added 
to the chiral fluctuations. To lowest order in density the 
QCD condensates give rise to a scalar self-energy Us = 
—(7NM/(m^f^)gs and a vector self-energy Ey — 4(m„ + 
md)M/(m%f%)Q- It is remarkable that the total self-energies, 
i.e. condensates plus chiral fluctuations, are very close to 
those obtained from DBHF calculations |12I23| . The re- 
sulting EOS is also shown in Fig. [Jin addition with that 
obtained after fine tuning to finite nuclei. Although the 
original EOS (case 1) is rather soft the inclusion of the 
condensates and the adjustment to finite nuclei results in 
an EOS with is finally stiff. 



2.3 Ab initio calculations 
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Fig. 2. Nuclear matter saturation points from relativistic 
(full symbols) and non-relativistic (open symbols) Brueckner- 
Hartree-Fock calculations based on different nucleon-nucleon 
forces. The diamonds show results from variational calcula- 
tions. Shaded symbols denote calculations which include 3- 
body forces. The shaded area is the empirical region of sat- 
uration. 



In ab initio calculations based on many-body tech- 
niques one derives the energy functional from first prin- 
ciples, i.e. treating short-range and many-body correla- 
tions cxplicitely. A typical example for a successful many- 
body approach is Brueckner theory |T2J. In the relativis- 
tic Brueckner approach the nucleon inside the medium 
is dressed by the self-energy S. The in-medium T-matrix 
which is obtained from the relativistic Bethe-Salpeter (BS) 
equation plays the role of an effective two-body interac- 
tion which contains all short-range and many-body cor- 
relations of the ladder approximation. Solving the BS- 
equation the Pauli principle is respected and intermediate 
scattering states are projected out of the Fermi sea. The 
summation of the T-matrix over the occupied states inside 
the Fermi sea yields finally the self-energy in Hartree-Fock 
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approximation. This coupled set of equations states a self- 
consistency problem which has to be solved by iteration. 

In contrast to relativistic DBHF calculations which 
came up in the late 80ties non-relativistic BHF theory 
has already almost half a century's history. The first nu- 
merical calculations for nuclear matter were carried out by 
Brueckner and Gammel in 1958 Despite strong efforts 
invested in the development of improved solution tech- 
niques for the Bethe-Goldstone (BG) equation, the non- 
relativistic counterpart of the BS equation, it turned out 
that, although such calculations were able to describe the 
nuclear saturation mechanism qualitatively, they failed 
quantitatively. Systematic studies for a large number of 
NN interactions were always allocated on a so-called Coester- 
line in the E/A—p plane which does not meet the empir- 
ical region of saturation. In particular modern one-boson- 
exchange (OBE) potentials lead to strong over-binding 
and too large saturation densities where relativistic cal- 
culations do a much better job. 

Fig. compares the saturation points of nuclear mat- 
ter obtained by relativistic Dirac-Brueckner-Hartree-Fock 
(DBHF) calculations using the Bonn potentials [SSj as 
bare NN interactions to non-relativistic Brueckner-Hartree- 
Fock calculations for various NN interactions. The DBHF 
results are taken from Ref. |23 (BM) and more recent 
calculations based on improved techniques are from |23j 
(Tubingen). Several reasons have been discussed in the 
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Fig. 3. Predictions for the EOS of symmetric nuclear matter 
from microscopic ab initio calculations, i.e. relativistic DBHF 
[23[ . non-relativistic BHF 18 and variational 15 calculations. 
For comparison also soft and hard Skyrme forces are shown. 



literature in order to explain the success of the relativistic 
treatment. The saturation mechanisms in relativistic and 
non-relativistic theories are quite different. In relativis- 
tic MFT the vector field grows linear with density while 
the scalar field saturates at large densities. The magni- 
tude and the density dependence of the scalar and vec- 
tor DBHF self-energy is similar to MFT, i.e. the single 
particle potential is the result of the cancellation of two 
large scalar and vector fields, each several hundred MeV 



in magnitude (see e.g. the effective mass in Fig. EJ. In 
BHF, on the other hand, the saturation mechanism takes 
place exclusively on the scale of the binding energy, i.e. a 
few ten MeV. It cannot be understood by the absence of a 
tensor force. In particular the second order 1-7T exchange 
potential (OPEP) is large and attractive at high densities 
and its interplay with Pauli-blocking leads finally to sat- 
uration. Relativistically the tensor force is quenched by 
a factor (m* D /M) 2 and less important for the saturation 
mechanism |37| . 

Three-body forces (3-BFs) have extensively been stud- 
ied within non-relativistic BHF [THj and variational calcu- 
lations JS] . The contributions from 3-BFs are in total re- 
pulsive which makes the EOS harder and non-relativistic 
calculations come close to their relativistic counterparts. 
The same effect is observed in variational calculations [TB] 
shown in Fig. The variational results shown contain 
boost corrections (Sv) which account for relativistic kine- 
matics and lead to additional repulsion JS]. Both, BHF 
[T%| and the variational calculations from ^1 are based 
on the latest AVis version of the Argonne potential. In 
both cases phenomenological 3-body-forces are used, the 
Tucson-Melbourne 3-BF in [TBI and the Urbana IX 3-BF 1 
in |15|. It is often argued that in non-relativistic treat- 
ments 3-BFs play in some sense an equivalent role as the 
dressing of the two-body interaction by in-medium spinors 
in Dirac phenomenology. Both mechanisms lead indeed 
to an effective density dependent two-body interaction V 
which is, however, of different origin. One class of 3-BFs 
involves virtual excitations of nucleon-antinucleon pairs. 
Such Z-graphs are in net repulsive and can be considered 
as a renormalization of the meson vertices and propaga- 
tors. A second class of 3-BFs is related to the inclusion 
of explicit resonance degrees of freedom. The most impor- 
tant resonance is the Z\(1232) isobar which provides at 
low and intermediate energies large part of the interme- 
diate range attraction. Intermediate A states appear in 
elastic NN scattering only in combination with at least 
two-isovector-meson exchange (7T7t, np, . . .). Such box di- 
agrams can satisfactorily be absorbed into an effective a- 
exchange |36|. The maintenence of explicit A DoFs gives 
rise to additional saturation, shifting the saturation point 
away from the empirical region |2l)|. However, as pointed 
out e.g. in |27| the inclusion of non-nucleonic DoFs has to 
be performed with caution: Freezing out resonance DoFs 
generates automatically a class of three-body forces which 
contains nucleon-resonance excitations. There exist strong 
cancellation effects between the repulsion due to box dia- 
grams and contributions from 3-BFs. Non-nucleonic DoFs 
and many-body forces should therefore be treated on the 
same footing. Such a treatment may be possible with the 
next generation of nucleon-nucleon forces based on chi- 
ral perturbation theory 



which allows a systematic 
generation of three-body forces. Next-to-leading order all 
3-BFs cancel while non-vanishing contributions appear at 
NNLO. 



1 Using boost corrections the repulsive contributions of the 
UIX interaction are reduced by about 40% compared to the 



original ones in [15) 
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Fig.|31compares the equations of state from the differ- 
ent approaches: DBHF from Ref. [S^] based the Bonn A in- 
teraction 2 36 , BHF 18 and variational calculations |15|. 
The latter ones are based on the Argonne AVis potential 
and include 3-body forces. All the approaches use modern 
high precision NN interactions and represent state of the 
art calculations. Two phenomenological Skyrme function- 
als which correspond to the limiting cases of a soft (K=200 
MeV) and a hard (K=380 MeV) EOS are shown as well. 
In contrast to the Skyrme interaction Q where the high 
density behavior is fixed by the compression modulus, in 
microscopic approaches the compression modulus is only 
loosely connected to the curvature at saturation density. 
DBHF Bonn A has e.g. a compressibility of K=230 MeV. 
Below 3po both are not too far from the soft Skyrme EOS. 
The same is true for BHF including 3-body forces. 

When many-body calculations are performed, one has 
to keep in mind that elastic NN scattering data constrain 
the interaction only up to about 400 MeV, which corre- 
sponds to the pion threshold. NN potentials differ essen- 
tially in the treatment of the short-range part. A model 
independent representation of the NN interaction can be 
obtained in EFT approaches where the unresolved short 
distance physics is replaced by simple contact terms. In 
the framework of chiral EFT the NN interaction has been 
computed up to N 3 LO 39 40 . An alternative approach 
which leads to similar results is based on renormalization 
group (RG) methods 0J- In the V"i ow k approach a low- 
momentum potential is derived from a given realistic NN 
potential by integrating out the high-momentum modes 
using RG methods. At a cutoff A ~ 2 fm _1 all the different 
NN potential models were found to collapse to a model- 
independent effective interaction V\ OV! k- When applied to 
the nuclear many-body problem low momentum interac- 
tions do not require a full resummation of the Brueck- 
ner ladder diagrams but can already be treated within 
second-order perturbation theory 05]. However, without 
repulsive three-body-forces isospin saturated nuclear mat- 
ter was found to collapse. Including 3-BFs first promising 
results have been obtained with VI ow k 021, however, nu- 
clear saturation is not yet described quantitativley. More- 
over, one has to keep in mind that - due to the high mo- 
mentum cut-offs - EFT is essentially suitable at moderate 
densities. 



3 EOS in symmetric and asymmetric nuclear 
matter 

Fig. El compares now the predictions for nuclear and neu- 
tron matter from microscopic many-body calculations - 
DBHF [2E1 and the 'best' variational calculation with 3- 
BFs and boost corrections [21 ~ to phenomenological ap- 
proaches and to EFT. As typical examples for relativistic 
functionals we take NL3 021 as one of the best RMF fits to 
the nuclear chart and a phenomenological density depen- 
dent RMF functional DD-TW from H2|. ChPT+corr. is 

2 The high density behavior of the EOS obtained with dif- 
ferent interaction, e.g. Bonn B or C is very similar. |23) 
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Fig. 4. EOS in nuclear matter and neutron matter. 
BHF/DBHF and variational calculations are compared to 
phenomenological density functionals NL3 and DD-TW and 
ChPT+corr.. The left panel zooms the low density range. 



based on chiral pion-nucleon dynamics including conden- 
sate fields and fine tuning to finite nuclei (case 3 in Fig. 
0. As expected the phenomenological functionals agree 
well at and below saturation density where they are con- 
strained by finite nuclei, but start to deviate substantially 
at supra-normal densities. In neutron matter the situa- 
tion is even worse since the isospin dependence of the 
phenomenological functionals is less constrained. The pre- 
dictive power of such density functionals at supra-normal 
densities is restricted. Ab initio calculations predict through- 
out a soft EOS in the density range relevant for heavy ion 
reactions at intermediate and low energies, i.e. up to about 
three times po- There seems to be no way to obtain an EOS 
as stiff as the hard Skyrme force shown in Fig.|3]or NL3. 
Since the nn scattering lenght is large, neutron matter at 
subnuclear densities is less model dependent. The micro- 
scopic calculations (BHF/DBHF, variational) agree well 
and results are consistent with 'exact' Quantum-Monte- 
Carlo calculations |29j. 




Fig. 5. Symmetry energy as a function of density as predicted 
by different models. The left panel shows the low density region 
while the right panel displays the high density range. 
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Y p which is 



In isospin asymmetric matter the binding energy is a 
functional of the proton and neutron densities, character- 
ized by the asymmetry parameter j3 = Y n 
the difference of the neutron and proton fraction Yi 
Pi/ p ,i — n,p. The isospin dependence of the energy func- 
tional can be expanded in terms of /3 which leads to a 
parabolic dependence on the asymmetry parameter 



E. 
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Fig. compares the symmetry energy predicted from the 
DBHF and variational calculations to that of the empir- 
ical density functionals already shown in Fig. 0] In ad- 
dition the relativistic DD-p5 RMF functional 01] is in- 
cluded. Two Skyrme functionals, SkM* and the more re- 
cent Skyrme-Lyon force SkLya represent non-relativistic 
models. The left panel zooms the low density region while 
the right panel shows the high density behavior of E aym . 
Remarkable is that most empirical models coincide around 



p ~ 0.6po where E s 



24 MeV. This demonstrates that the non-relativistic effective mass ; 



functional varying thereby the linear asymmetry param- 
eter a 4 from 30 to 38 MeV |55|. In it was concluded 
that charge radii, in particular the skin thickness r n — r p in 
heavy nuclei constrains the allowed range of a 4 to 32 -j- 36 
MeV for relativistic functionals. 

Fig. H3 displays the correlation between the skin thick- 
ness in 208 Pb and a 4 obtained within various models. The 
skin thickness depends, however, not only on the symme- 
try energy but there exists a close correlation between a 4 
and the compression modulus K |35| . This correlation is 
of importance when these quantities are extracted from 
finite nuclei (see discussion by Shlomo et al. in Sec. 1.4). 



3.1 Effective nucleon masses 

The introduction of an effective mass is a common concept 
to characterize the quasi-particle properties of a particle 
inside a strongly interacting medium. In nuclear physics 
exist different definitions of the effective nucleon mass 
which are often compared and sometimes even mixed up: 

and the relativistic 



constraints from finite nuclei are active for an average 
density slightly above half saturation density. However, 
the extrapolations to supra-normal densities diverge dra- 
matically. This is crucial since the high density behavior 
of -E sym is essential for the structure and the stability of 
neutron stars (see also the discussion in Sec. V.5). The 
microscopic models show a density dependence which can 
still be considered as asy-stiff. DBHF 26 is thereby stiffcr 
than the variational results of (T^l- The density depen- 
dence is generally more complex than in RMF theory, in 
particular at high densities where E sym shows a non-linear 
and more pronounced increase. Fig.[S]clearly demonstrates 
the necessity to constrain the symmetry energy at supra- 
normal densities with the help of heavy ion reactions. The 
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Dirac mass m* D . These two definitions are based on dif- 
ferent physical concepts. The nonrelativistic mass param- 
eterizes the momentum dependence of the single-particle 
potential. The relativistic Dirac mass is defined through 
the scalar part of the nucleon self-energy in the Dirac 
field equation which is absorbed into the effective mass 
m* D = M + Ss(k, Uf)- The Dirac mass is a smooth func- 
tion of the momentum. In contrast, the nonrelativistic ef- 
fective mass - as a model independent result - shows a 
narrow enhancement near the Fermi surface due to an 
enhanced level density 05]. For a recent review on this 
subject and experimental constraints on m* NR see |46| . 

While the Dirac mass is a genuine relativistic quan- 
tity the effective mass rn* NR is determined by the single- 
particle energy 
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is a measure of the non-locality of the single-particle 
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Fig. 6. Skin thickness in 208 Pb versus the linear symmetry 
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hatched area in Fig.[5]displays the range of -E S ym which has 
been obtained by constructing a density dependent RMF 
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Fig. 7. The effective mass in isospin symmetric nuclear matter 
as a function of the momentum k at different densities deter- 
mined from relativistic Brueckner calculations. 
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potential U (real part) which can be due to non-localities 
in space, resulting in a momentum dependence, or in time, 
resulting in an energy dependence. In order to clearly 
separate both effects, one has to distinguish further be- 
tween the so-called k-mass and the E-mass ri7J. The spa- 
tial non-localities of U are mainly generated by exchange 
Fock terms and the resulting k-mass is a smooth func- 
tion of the momentum. Non-localities in time are gener- 
ated by Brueckner ladder correlations due to the scatter- 
ing to intermediate states which are off-shell. These are 
mainly short-range correlations which generate a strong 
momentum dependence with a characteristic enhancement 
of the E-mass slightly above the Fermi surface [45117147] . 
The effective mass defined by Eq. (jSJ) contains both, non- 
localities in space and time and is given by the product 
of k-mass and E-mass JT7J. In Fig. [7J the nonrelativis- 
tic effective mass and the Dirac mass, both determined 
from DBHF calculations ^Bli are shown as a function 
of momentum k at different Fermi momenta of kp = 
1.07, 1.35, 1.7 fm . 

m NB. snows the typical peak struc- 
ture as a function of momentum around kp which is also 
seen in BHF calculations 07] • The peak reflects the in- 
crease of the level density due to the vanishing imaginary 
part of the optical potential at kp which is also seen, e.g., 
in shell model calculations |45I17| . One has, however, to 
account for correlations beyond mean field or Hartree- 
Fock in order to reproduce this behavior. Fig. |H| com- 



0.2 




_ 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 _ 
: nonrelativistic mass : 


_ i i i i i i i 

- Dirac mass 

- *' 


i 

— DBHF (Bonn A) 
□ BHF (Bonn A) 

— QHD-I 
■-■ NL3 
■■■■ DD-TW 




- S%v - 



2 3 1 

P / P n P / P, 



Fig. 8. Nonrelativistic and Dirac effective mass in isospin sym- 
metric nuclear matter as a function of the density for various 
models. 



pares the density dependence of the two effective masses 
determined at kp. Both masses decrease with increasing 
density, the Dirac mass continously, while rn* NR starts to 
rise again at higher densities. Phcnomenological density 
functionals (QHD-I, NL3, DD-TW) yield systematically 
smaller values of m* NR than the microscopic approaches. 
This reflects the lack of nonlocal contributions from short- 
range and many-body correlations in the mean field ap- 
proaches. 



3.1.1 Proton-neutron mass splitting 



A heavily discussed topic is in the moment the proton- 
neutron mass splitting in isospin asymmetric nuclear mat- 
ter. This question is of importance for the forthcoming 
new generation of radioactive beam facilities which are de- 
voted to the investigation of the isospin dependence of the 
nuclear forces at its extremes. However, presently the pre- 
dictions for the isospin dependences differ substantially. 
BHF calculations |18I47| predict a proton-neutron mass 
splitting of m* NRn > m* NRp . This stands in contrast to 
relativistic mean-field (RMF) theory. When only a vector 
isovector p-meson is included Dirac phenomenology pre- 
dicts equal masses m* D n = m* D while the inclusion of the 
scalar isovector 5-meson, i.e. p + S, leads to m* D n < m* D 
|44| . When the effective mass is derived from RMF the- 
ory, it shows the same behavior as the corresponding Dirac 
mass, namely m* NR n < m* NRp 44 j. Conventional Skyrme 
forces, e.g. SkM*, lead to m* NRn < m* NRp gUj while 
the more recent Skyrme-Lyon interactions (SkLya) pre- 
dict the same mass splitting as RMF theory. The pre- 
dictions from relativistic DBHF calculations are in the 
literature still controversial. They depend strongly on ap- 
proximation schemes and techniques used to determine 
the Lorentz and the isovector structure of the nucleon self- 
energy. In the approach originally proposed by Brockmann 
and Machleidt |U one extracts the scalar and vector self- 
energy components directly from the single-particle po- 
tential. Thus, by a fit to the single-particle potential mean 
values for the self-energy components are obtained where 
the explicit momentum-dependence has already been aver- 
aged out. In symmetric nuclear matter this method is rel- 
atively reliable but the extrapolation to asymmetric mat- 
ter is ambiguous [21] ■ Calculations based on this method 
predict a mass splitting of m* D n > m* D 50 . On the 
other hand, the components of the self-energies can di- 
rectly be determined from the projection onto Lorentz 
invariant amplitudes 20 22 23 24 26 51 . Projection tech- 
niques are involved but more accurate and yield the same 
mass splitting as found in RMF theory when the S -meson 
is included, i.e. m* D n < m* D 22 2 4l26j . Recently also the 
non-relativistic effective mass has been determined with 
the DBHF approach and here a reversed proton-neutron 
mass splitting was found, i.e. m* NRn > m* NRp 0Bj- Thus 
DBHF is in agreement with the results from nonrelativis- 
tic BHF calculations. 

Experimentally accessable is thep— n mass splitting, or 
the magnitude of the corresonding isovector effective mass 
m v , {Sr = M 1 - ~ 7^ — ) through the electric dipole 

V NR NR,n 

photoabsorption cross section, i.e. through an enhance- 
ment of the Thomas-Reiche-Kuhn sum rule by the factor 
m/rriy. However, values derived from GDR measurements 
range presently from m v /m = 0.7 -j- 1.05 46 52I53 | . The 
forthcoming radioactive beam facilitites will certainly im- 
prove on this not yet satisfying situation. 
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3.2 Optical potentials 

The second important quantity related to the momen- 
tum dependence of the mean field is the optical nucleon- 
nucleus potential. At subnormal densities the optical po- 
tential £/opt is constraint by proton-nucleus scattering data 
|54| and at supra-normal densities constraints can be de- 
rived from heavy ion reactions {5M5BE7j . In a relativis- 
tic framework the optical Schroedinger-equivalent nucleon 
potential (real part) is defined as 



U op t — —£s + jjJ-^V 



r 2 - r 2 

2M 



(6) 



One should thereby note that in the literature sometimes 
also an optical potential, given by the difference of the 
single-particle energies in medium and free space U = 
E — V M 2 + k 2 is used jSS] which should be not mixed 
up with In a relativistic framework momentum inde- 
pendent fields Ss,V ( as e -g- i n RMF theory) lead always 
to a linear energy dependence of £/ pt- As seen from Fig. El 
DBHF reproduces the empirical optical potential [HI] ex- 
tracted from proton-nucleus scattering for nuclear matter 
at po reasonably well up to a laboratory energy of about 
0.6-0.8 GeV. However, the saturating behavior at large 
momenta cannot be reproduced by this calculations be- 
cause of missing inelasticities, i.e. the excitation of isobar 
resonances above the pion threshold. When such contin- 
uum excitations are accounted for optical model cacula- 
tions are able to describe nucleon-nucleus scattering data 
also at higher nergies |58j . In heavy ion reactions at in- 
cident energies above 1 AGeV such a saturating behav- 
ior is required in order to reproduce transverse flow ob- 
servables [57]. One has then to rely on phenomenological 
approaches where the strength of the vector potential is 
artificially suppressed, e.g. by the introduction of addi- 
tional form factors |57) or by energy dependent terms in 
the QHD Lagrangian [5J] (D 3 C model in Fig© . 

The isospin dependence, expressed by the isovector 
optical potential L7 iso = (U op t, n — U op t, p ) / (2(3) is much 
less constrained by data. The knowledge of this quan- 
tity is, however, of high importance for the forthcoming 
radioactive beam experiments. The right panel of Fig. 
E] compares the predictions from DBHF 26 and BHF 
|6D] to the phenomenological Gogny and Skyrme (SkM* 
and SkLya) forces and a relativistic T — p approximation 
[H2] based on empirical NN scattering amplitudes [SH] ■ At 
large momenta DBHF agrees with the tree-level results 
of While the dependence of U- mo on the asymmetry 
parameter (5 is found to be rather weak |26I60| . the pre- 



dicted energy and density dependences are quite differ- 
ent, in particular between the microscopic and the phe- 
nomenological approaches. The energy dependence of i7i so 
is very little constrained by data. The old analysis of op- 
tical potentials of scattering on charge asymmetric tar- 
gets by Lane |M] is consistent with a decreasing potential 
as predicted by DBHF/BHF, while more recent analyses 
based on Dirac phenomenology [53] come to the opposite 
conclusions. RMF models show a linearly increasing en- 
ergy dependence of L7 iso (i.e. quadratic in k) like SkLya, 
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Fig. 9. Nucleon optical potential in nuclear matter at po- On 
the left side DBHF calculations for symmetric nuclear mat- 
ter from [20] and [23 

are compared to the phenomenological 
models NL3 and D 3 C |59| and to the p-A scattering analysis 
of |54) . The right panel compares the iso- vector optical poten- 
tial from DBHF |55] and BHF [§U] to phenomenological RMF 
|61) , Gogny and Skyrme forces and to a relativistic T — p 
approximation |fi2| . 



however generally with a smaller slope (see discussion in 
|44j ) . To clarify this question certainly more experimental 
efforts are necessary. 



4 Transport models 

The difficulty to extract information on the EOS from 
heavy ion reactions lies in the fact that the colliding sys- 
tem is over a large time span of the reaction out of global 
and even local equilibrium. At intermediate energies the 
relaxation time needed to equilibrate coincides more or 
less with the high density phase of the reaction. Hence, 
non-equilibrium effects are present all over the compres- 
sion phase where one essentially intends to study the EOS 
at supra- normal densities. Experimental evidences for in- 
complete equilibration even in central collisions have been 
found by isospin tracing of projectile and target nuclei |66| 
and by different variances of longitudinal and transverse 
rapidity distributions }67j. To account for the temporal 
space time evolution of the reactions requires dynamical 
approaches which are based on kinetic transport theory. 
In the following we briefly discuss the various approaches 
which are mainly used in order to describe the reaction 
dynamics at low and intermediate energies. 



4.1 Boltzmann-type kinetic equations 

The theoretical basis for the description of the collision dy- 
namics at energies ranging from the Fermi regime up to 
1-2 AGeV is hadronic non-equilibrium quantum transport 
field theory |68j . The starting point of non-equilibrium 
QFT is the Schwinger-Keldysh formalism for many-body 
Green functions in non-equilibrium configurations. The 
one-body Green function is defined as the expectation 
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value of the time ordered product of fermionic field op- 
erators G(1,I') = H) < T sk (&(l)&(l')) > where T sk 
defines the temporal sequence of the field operators. In 
non-equilibrium time reversal invariance is violated and 



thus the application of T s k 
nations 

G c = -i(T c [V(l)W(l')}), 



which describes the phase space evolution of the 1-particle 
distribution /(x, k, t) under the influence of the mean field 
(which enters via the real part of the self energy, i.e. 
via m* D = M - S s and k*^ = - S^) and binary 



leads to four possible combi- collisions determined by the transition amplitude W 



G a - 
G< 



i{WW)) (7) 



where T c (T°) is the causal (anti-causal) time ordering op- 
erator. The physical quantity of interest is the correlation 
function G < since it corresponds in the equal time limit 
to the density lim ti ,_ tl G < (1, 1') = (+i)p(x\, x^, f). How- 
ever, the four Green-functions are related through equa- 
tions of motion (Kadanoff-Baym Equations) for the corre- 
lation G <>> and the retarded and advanced G^ functions 
(the retarded and advanced Green functions are defined 
via G+- = G c -G<<> = G > - < -G a ). From the Kadanoff- 
Baym Eqs. one obtains a kinetic equation for the correla- 
tion function G < 



m^\T{hk2\kski)\ ■ Final state Pauli blocking is accounted 
for by the blocking factors (1—fA in with fi = /(x, kj, i). 
The physical parameters entering into the kinetic equation 
are the mean field, i.e. the nuclear EOS, and elementary 
cross sections for 2-particle scattering processes. Thus one 
can test the high density behavior of the nuclear EOS 
in heavy ion collisions and the in-medium modifications 
of cross sections, which also influence the stopping prop- 
erties of the colliding system. Above the pion threshold 
where inelastic processes start to play an important role 
the Eq. @ becomes a coupled channel problem for nucle- 
onic, nucleon resonance and mesonic degrees of freedom. 
The collision integral, i.e. the right hand side of Eq. © 
has to be extended for the corresponding inelastic and ab- 
sorptive processes and the new degrees of freedom must be 
propagated in their mean fields. In practice the transport 



DG< - G < D* - (ReS + G < 



-ReG + E<' 



= - 2 (s'c,< 



D = —id Xl /2M is the Schrodinger operator or in a rela- 
tivistic framework the Dirac operator (D = ij^d Xl — M) 
and E < ' > ' ± are the self energies. The introduction of re- 
tarded and advanced functions allows to interpret the real 
part of the retarded self energy as a mean field while the 
imaginary part describes the absorption or finite life times 
of quasi particles (dressed nucleons) [HHj- The self-energy 
E contains all higher order correlations and couples the 
one-body kinetic equation (jHJ) to the corresponding equa- 
tions for the two- and 3-body densities and so forth. This 
requires to truncate the Dyson-Schwinger hierarchy which 
is usually done at the two-body level and leads to the lad- 
der approximation for T-matrix, i.e. the Bethe-Salpeter 
equation. 

The formal structure of the kinetic equation JHJ is 
complex and one should solve JHJ) together with the cor- 
responding kinetic equations for G ± which describe the 
spectral properties of the phase space distribution. Simul- 
taneously the self-energies should be derived for arbritrary 
non-equilibrium situations (HHj • A solution of the full self- 
consistency problem has not yet achieved. In praxis one 
applies further approximations. The most important ones 
are the gradient expansion (a semi-classical approxima- 
tion to first order in K) and the quasi-particle approxima- 
tion which sets the particles on mass shell. The result is a 
Boltzmann-type transport equation, which is known as the 
Boltzmann-Uheling-Uhlcnbcck (BUU) transport equation 
|69|. In its relativistic form the (R)BUU equation reads 



< _|_ equation is solved within the testparticle method which 

— G ReE ) — (E ReG describes the phase space distribution / as an incoherent 
q<jj> _ £<q> _ q> jj<\ C8"f um °f P om tTike quasi particles [HI] or static gaussians 
^ ^ ' * -]7D| which propagate on classical trajectories. Relativistic 

formulations of the two methods were developed in |7I| 
and [El- 



4.2 Quantum Molecular Dynamics (QMD) 

An alternative approach to the kinetic BUU equation is 
Quantum Molecular Dynamics (QMD) |73I74I75| . QMD 
is a N-body approach which simulates heavy ion reactions 
on an event by event basis taking fluctuations and corre- 
lations into account. The QMD equations are formally de- 
rived from the assumption that the N-body wave function 
<P can be represented as the direct product of single co- 
herent states <P — Yli <l>i which are described by Gaussian 
wave packets. Anti-symmetrization is not taken into ac- 
count. A Wigner transformation yields the corresponding 
phase space representation of The equations of motion 
of the many-body system are obtained by the variational 
principle starting from the action S = J #*] (with 
the Lagrangian functional C =< ^\ih^ — H\<P >). The 
Hamiltonian H contains a kinetic contribution and mu- 
tual two-body interactions Vij. The variational principle 
leads finally to classical equations of motion for the gener- 
alized coordinates and k^ of the Gaussian wave packets 



E 



< Vu >= V k , < H > 



'D 



-k™dZkl)dl-{m* D d»m,* D 
d 3 k 2 d 3 k 3 d 3 k 4 



k* u d%K 



< Va >= V aj < H > 



El Et 



El 

K4 



W S 4 {k + k 2 -k 3 - fc 4 ) 



2(2tt) 8 

hh (i - /) (i - / 2 ) - fh (i - h) (i - U) 



The two-body interaction Vij can be e.g. taken from BHF 
calculations |75| or from local Skyrme forces which are 
usually supplemented by an empirical momentum depen- 
(9) dence in order to account for the energy dependence of the 
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optical nucleoli- nucleus potential \73\ . Binary collisions 
are treated in the same way as in BUU models. Further- 
more there exist relativistic extensions, i.e. RQMD and 
the UrQMD model which has been developed to simulate 
heavy ion collisions at ultra-relativistic energies |74I76| . 

4.3 Antisymmetrized Molecular Dynamics 
(AMD/FMD) 

An extension of QMD, in particular designed for low ener- 
gies, are the Antisymmetrized Molecular Dynamics (AMD) 
|77) and Fermionic Molecular Dynamics (FMD) approaches 
|78j . In contrast to conventional QMD, the interacting 
system is represented by an antisymmetrized many-body 
wave function consisting of single-particle states which are 
localized in phase space. The equations of motion for the 
parameters characterizing the many-body state (e.g. po- 
sition, momentum, width and spin of the particles) are 
derived from a quantum variational principle. The models 
are designed to describe ground state properties of nuclei 
as well as heavy ion reactions at low energies (see also 
Chap. 1.2). 

4.4 Off-shell transport 

Essential for the validity of the classical equations of mo- 
tion is the quasi-particle approximation (QPA) which as- 
sumes that the spectral strength of a hadron is concen- 
trated around its quasi-particle pole. Particle widths can, 
however, dramatically change in a dense hadronic envi- 
ronment. To first order in density the in- medium width 
of a hadron in nuclear matter can be estimated by the 
collision width r tot = T vac + r co11 , r coU = jvap B with 
v the hadron velocity relative to the surrounding matter 
and a the total hadron-nucleon cross section. A consistent 
treatment of the off-shell dynamics, i.e. a solution of the 
quantum evolution equations for the correlation functions 
G <,> has up to now only been performed for toy models 
and simplified geometries |79I8()| or in first order gradi- 
ent approximation leading to an extended quasi-particle 
picture [SJ. Comparing the nonlocal extension of BUU 
with standard simulations a visible effect of nonlocal cor- 
relations is seen and a better agreement with measured 
charge density distributiond jH2 or particle spectra 
due to the virial corrections has been found. To develop a 
consistent lattice quantum transport for non-uniform sys- 
tems and realistic interactions will be one of the future 
challenges in theoretical heavy ion physics. 

On the other hand, substantial progress has been made 
in the recent years to map part of the off-shell dynam- 
ics on a modified test-particle formalism |84I85| . This al- 
lows to apply off-shell dynamics, although in a simplified 
form, to the complex space time evolution of a heavy ion 
reaction. The present knowledge of off-shell matrix ele- 
ments is, however, rather limited and theoretical investi- 
gations are scarce [EE]. The off-shell T- matrix has been 
used in order to calculate the duration and non-locality 
of a nucleon-nucleon collision [HZj- The question to what 



degree a depletion of the Fermi surface due to particle- 
hole excitations and the high momentum tails of the nu- 
clear spectral functions will affect subthreshold particle 
production is not so obvious to answer. The high mo- 
mentum tails correspond to deeply bound states which 
are off-shell and to treat such states in a standard trans- 
port approach like on-shell quasi-particles would violate 
energy-momentum conservation. Energy-momentum con- 
servation can be achieved consistently by the nonlocal ki- 
netic theory 88 taking into account first order off-shell 
effects. The contribution of the nuclear short-range corre- 
lations to subthreshold K + production in p+A reactions 
have e.g. been estimated in 89 . The removal energy for a 
high momentum state compensates the naively expected 
energy gain and the short-range correlations do therefore 
not significantly contribute to subthreshold particle pro- 
duction [55]. The situation changes, however, when the 
medium is heated up and high momentum particles be- 
come on-shell or when the spectral distributions of the 
produced hadrons themselves are broadened. 



5 Constraints from heavy ion collisions 
5.1 Flow and stopping 

On of the most important observable to constrain the nu- 
clear forces and the underlying EOS at supra-normal den- 
sities is the collective nucleon flow |5U]. It can be charac- 
terized in terms of anisotropies of the azimuthal emission 
pattern. Expressed in terms of a Fourier series 

dN 

— cx 1 + 2v 1 cos((/)) + 2v 2 cos(2(/)) + ... (10) 
dip 

this allows a transparent interpretation of the coefficients 
V\ and v 2 . The dipole term v\ arises from a collective side- 
ward deflection of the particles in the reaction plane and 
characterizes the transverse flow in the reaction plane. The 
second harmonics describes the emission pattern perpen- 
dicular to the reaction plane. For negative v 2 one has a 
preferential out-of-plane emission. The phenomenon of an 
out-of-planc enhancement of particle emission at midra- 
pidity is called squeeze-out. 

The transverse flow v\ has been found to be sensitive to 
the EOS and, in particular in peripheral reactions, to the 
momentum dependence of the mean field {25EHI- The el- 
liptic flow v 2 , in contrast, is very sensitive to the maximal 
compression reached in the early phase of a heavy ion reac- 
tion. The cross over from preferential in-plane flow v 2 < 
to preferential out-off-plane flow v 2 > around 4-6 AGeV 
has also led to speculations about a phase transition in 
this energy region which goes along with a softening of the 
EOS The present situation between theory and exper- 
iment is illustrated in Fig. EH (from The BUU studies 
from Danielewicz et al. and the Giessen group (Larionov 
et al.) investigated the EOS dependence while Persram 
et al. find a sensitivity of v 2 on the medium dependence 
of the NN cross sections. Finally, non-equilibrium effects 
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Fig. 10. Elliptic flow excitation function at SIS energies. Var- 
ious theoretical studies using different EOS's ((a),(b)), or dif- 
ferent cross sections (c) or DBHF mean fields in the LDA ap- 
proach and including further non-equilibrium effects (ACNM) 
are compared to FOPI data (symbols). Figure is taken from 



have been investigated at the level of the effective inter- 
action in |56I93| . It has been found that the local phase 
space anisotropies of the pre-equilibrium stages of the re- 
actions reduce the repulsion of the mean field and soften 
the corresponding EOS which allows a good description 
of the V2 data using microscopic DBHF mean fields (Gai- 
tanos et al.). However, Fig. ^] also demonstrates that v-i 
is generated by the interplay of the mean field and binary 
collisions which makes it difficult to extract exclusive in- 
formation on the EOS from the data. Here certainly fur- 
ther going studies are required. 

The next figure lTTl based on the studies of Danielewicz 
|94| , summarizes the status obtained within this model in 
terms of a band that represents the constraints from col- 
lective flow data. It based on a compilation from the anal- 
ysis of sideward and elliptic anisotropies, ranging from 
low SIS (Si a b ^ 0.2 2 AGeV) up to top AGS energies 
(Slab — 2 -=-11 AGeV), has been studied in [<2]. The con- 
clusion from this study was that, both, super-soft equa- 
tions of state (K=167 MeV) as well as hard EOSs (K>300 
MeV) are ruled out by data. At SIS energies existing flow 
data are consistent with a soft EOS [95155) . e.g. the The 
soft Skyrme EOS. In the models used by Danielewicz et 
al. |55I94| sideward flow favors a rather soft EOS with 
K=210 MeV while the development of the elliptic flow 
requires slightly higher pressures. The BHF and varia- 
tional calculations including 3-body-forces 3 fit well into 



the constrained area up to 4pn. At higher densities the 
microscopic EOSs, also DBHF, tend to be too repulsive. 
However, conclusions from flow data are generally com- 
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For the BHF + 3-BF calculation the pressure shown in Fig. 
Illl has been determined from the parameterization given in |19| 



Fig. 11. Constraints on the nuclear EOS from heavy ion flow 
data. The shaded area shows the pressure-density which is com- 
patible with heavy ion flow data according the analysis on |t)4) . 
The equations-of-state from the models shown in Fig.|3]are dis- 
played. 



plicated by the interplay of the compressional part of the 
nuclear EOS and the momentum dependence of the nu- 
clear forces. A detailed comparison to v\ and v-i data from 
FOPI [5fif9T] and KaoS [HH| for v x and v 2 below 1 AGeV 
favors again a relatively soft EOS with a momentum de- 
pendence close to that obtained from microscopic DBHF 
calculations 55 56 99 . In Fig. El the microscopic DBHF 
EOS (K=230 MeV) lies at the upper edge of the bound- 
ary, but is still consistent in the density range tested at SIS 
energies, i.e. up to maximally 3 po- This fact is further con- 
sistent with the findings of Gaitanos et al. 56 9(J| where 
a good description of v% and vi data at energies between 
0.2 and 0.8 AGeV has been found in RBUU calculations 
based on DBHF mean fields. As pointed out in [56193199] 
it is thereby essential to account for non-equilibrium ef- 
fects and the momentum dependence of the forces which 
softens the EOS compared to the equilibrium case (shown 
in Fig. Ell . 

As can be seen from Fig. llUI not only the nuclear EOS, 
but also the cross sections for elementary 2-particle scat- 
tering influences the collective dynamics, in particular the 
degree of stopping and hence the maximal compression 
achieved in the fireball region. A challenge is in this con- 
text the quantitative understanding of the recently ob- 
served strong correlations between maximal side flow v\ 
and maximal stopping in both excitation functions |67| 
(see Chap. 1.5 by W. Reisdorf et al.). Most collective flow 
analyses performed so far have were based on free cross 
sections, which works astonishingly well from a practi- 
cal point of view. However, within a consistent picture 

which is based on the Urbana IX 3-BF different to that used 
in ES]. 



12 



C. Fuchs, H.H. Wolter: Modelization of the EOS 



one should treat the in-medium effects in both, real (nu- 
clear EOS) and the imaginary part (cross sections) of 
the interaction on the same footing. Many-body calcula- 
tions (BHF/DBHF) predict an essential reduction of the 
elastic NN cross section with increasing baryon density 
[751100185] and in |101| a similar reduction was proposed 
for the inelastic channels in oder to describe pion multi- 
plicities in the 1-2 AGeV region. One can therefore expect 
observable signals in heavy ion collisions. In fact, recent 
QMD studies of stopping and transparency observables 
have shown that the data can be reproduced when the 
free cross section is reduced by a factor of 0.5 [102] . These 
findings are supported by transport calculations using mi- 
croscopic in- medium cross sections [1031104] . Therefore, 
for a reliable extraction of the high density nuclear EOS 
one should account for in-medium effects not only in the 
potential, but also in the cross sections. 

5.1.1 Isospin dependence of the EOS 

Another important aspect of heavy ion collisions is the 
investigation of the density dependence of the EOS for 
asymmetric matter. There exist abundent studies on this 
sector, either non-relativistically or relativistically. 

The momentum dependence of the isovector potential, 
Fig.lHl which is also closely related to the proton/neutron 
mass splitting of both, the non-relativistic rn* NR and the 
Dirac m* D effective mass, is one of the key questions which 
can be addressed by nuclear reactions induced by neutron- 
rich nuclei at RIA energies. Transverse and elliptic flow 
pattern as well p/n rapidity dsitributions have been sug- 
gested as possible observables to investigate the momen- 
tum dependence and the p/n mass splitting [10511061107] . 

Promising observables to pin down the density depen- 
dence of the symmetry energy are the iso-scaling behavior 
of fragment yields and the isospin diffusion in asymmetric 
colliding systems. In both cases recent NSCL-MSU data 
in combination with transport calculations are consistent 
with a value of E sym « 31 at po and rule out extremely " 
stiff " and " soft " density dependences of the symmetry 
energy |108I109| (see also Chap. III. 2). The same value 
has been extracted |110| from low energy elastic and (p,n) 
charge exchange reactions on isobaric analog states, i.e. 
p( 6 iJe, 6 Li*)n measured at the HMI. Such a behavior is 
also consistent with the predictions from many-body the- 
ory 15 2BJ- Also the p/n ratio at mid-rapidity has been 
found to be sensitive to the high density behavior of nu- 
clear symmetry energy |111| . 

In relativisitic approaches large attractive scalar and 
repulsive vector fields are required by Dirac phenomenol- 
ogy in order to describe simultaneously the central poten- 
tial and the strong spin-orbit force in finite nuclei [819112] . 
The situation is, however, less clear in the iso-vector sec- 
tor. There exists different possibilities to reproduce the 
same value of the 04 coefficient @: (a) by only an iso- 
vector vector p field like in most RMF models (NL3 etc.), 
or (b) by accounting for an additional iso-vector scalar 8 
field. Due to competing effects between attractive (scalar 
8) and repulsive (vector p) both alternatives can be fitted 



to the same empirical 04 parameter, but the inclusion of 
8 field leads to an essentially different high density be- 
havior of the symmetry energy [112] , The scalar 8 field 
is suppressed at high densities, whereas the vector field 
is proportional to the baryon density which makes the 
symmetry stiffer energy at supra-normal densities. Recent 
transport studies have shown that these subtle relativistic 
effects can be observed in the intermediate energy range by 
means of collective isospin flow, particle ratios and imbal- 
ance ratios of different particle species (protons, neutrons, 
pions and kaons) |112llU7fl3| . However, due to the lack of 
precise experimental data no definitive conclusions could 
be made so far. 



5.2 Particle production 

5.2.1 Pions 

With the start of the first relativistic heavy ion programs 
the hope was that particle production would provide a 
direct experimental access to nuclear EOS 113. At two 
times saturation density which is reached in the partici- 
pant zone of the reactions without additional compression 
the difference between the soft and hard EOS shown in 
Fig. is about 13 MeV in binding energy. If the matter is 
compressed up to 3po the difference is already ~ 55 MeV. 
It was expected that the compressional energy should be 
released into the creation of new particles, primarily pi- 
ons, when the matter expands [113] . However, pions have 
large absorption cross sections and they turned out not 
to be suitable messengers of the compression phase. They 
undergo several absorption cycles through nucleon reso- 
nances (Ntt <-> A) and freeze out at final stages of the 
reaction and at low densities. Hence pions loose most of 
their knowledge on the compression phase and are not 
really sensitive probes for stiffness of the EOS [114] , How- 
ever, they carry information on the isotopic composition 
of the matter which is to some extent conserved until 
freeze-out. The final tt~ /n + ratio was found to be sen- 
sitive to the initial n/p composition of the matter which, 
on the other hand, is influenced by the isospin dependence 
of the nuclear forces [1151116) . In O] a reduction of the 
7r~/7r + -ratio was found when the J-meson was included in 
the RMF approach. The effects are, however, moderate, 
i.e. at the 10-20% level, and most pronounced at extreme 
phase space regions, e.g. at the high energy tails of pt 
spectra 116 I61I117| . Systematic measurements, e.g. from 
the FOPI Collaboration may help to constrain the isospin 
dependence by pionic observables. 

5.2.2 Kaons 

After pions turned out to fail as suitable messengers, K + 
mesons were suggested as promising tools to probe the nu- 
clear EOS, almost 20 years ago [119] , The cheapest way 
to produce a K + meson is the reactions NN — ► NAK + 
which has a threshold of £q a b = 1.58 GeV kinetic energy 
for the incident nucleon. When the incident energy per 
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Fig. 12. K + multiplicities in inclusive C+C, Ni+Ni, Au+Au 
and C+Au reactions at 1 AGeV. QMD calculations using a 
hard/soft nuclear EOS are compared to KaoS data |118| . Fig- 
ure is taken from II 1? 
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nuclcon in a heavy ion reactions is below this values one 
speaks about subthreshold kaon production. Subthreshold 
kaon production is in particular interesting since it ensures 
that the kaons originate from the high density phase of the 
reaction. The missing energy has to be provided either by 
the Fermi motion of the nucleons or by energy accumu- 
lating multi-step reactions. Both processes exclude signifi- 
cant distortions from surface effects if one goes sufficiently 
far below threshold. In combination with the long mean 
free path subthreshold K + production is an ideal tool to 
probe compressed nuclear matter in relativistic heavy ion 
reactions. 

Already in the first theoretical investigations by trans- 
port models it was noticed the K + yield reacts rather sen- 
sitive on the EOS |12UI121I122] . Both, in non-relat ivistic 
QMD calculatio ns based on soft/hard Skyrme forces |12Ull2lilZj] 4 \- 
and in RBUU |122I124| with soft/hard versions of the 
(non-linear) crcj-model the K + yield was found to be about 
a factor 2-3 larger when a soft EOS is applied compared 
to a hard EOS. At that time the available data favored a 
soft equation of state |121I122I124| . However, at that stage 
the theoretical calculations were still burdened with large 
uncertainties. First of all, it was noticed 120 121 that 
the influence of the repulsive momentum dependent part 
of the nuclear interaction leads to a strong suppression of 
the kaon abundances which made a quantitative descrip- 
tion of the available data more difficult. Moreover, at that 
time the pion induced reaction channels ttB — ► YK + 
have not yet been taken into account. These additional 
channels which contribute up to 30 + 50% to the total 
yield enabled to explain the measured yields with realis- 
tic momentum dependent interactions |123I125| . A break- 
through was achieved when the COSY-11 Collaboration 
measured the pp — > pK + A reactions at threshold |126j 
which constrains the strangeness production cross sections 
NN — ► NK+Y Within the last decade the KaoS Col- 
laboration has performed systematic measurements of the 
K + production far below threshold 1 1411271128) . Based 
on the new data situation, the question if valuable infor- 
mation on the nuclear EOS can be extracted has been 
revisited and it has been shown that subthreshold K + 



production provides indeed a suitable and reliable tool for 
this purpose |12^ll3()ll3l| . 

Fig. ^] compares measured the K + multiplicities as a 
function of participating nucleons A valt in Au+Au, Ni+Ni, 
C+Au and C+C reactions at 1 AGeV to QMD calcu- 
lation using a soft/hard momentum dependent Skyrme 
force |118| . This figure demonstrates thereby the inter- 
play between ^4 par t, system size and EOS. A significant 
dependence of the kaon multiplicities on the nuclear EOS 
requires a large amount of collectivity which is easiest 
reached in central reactions of heavy mass systems. Conse- 
quently, the EOS dependence is most pronounced in cen- 
tral Au+Au reactions. Also in Ni+Ni effects are still siz- 
able while the small C+C system is completely insensitive 
on the nuclear EOS even in most central reactions. The 
data available for Au+Au and Ni+Ni support the soft 
EOS. Interesting is in this context the asymmetric C+Au 
system: Though in central C+Au reactions the number of 
participants is comparable to Ni+Ni the K + yield does 
not depend on the EOS. This indicates again that a sen- 
sitivity on the EOS is not only a question of A paTt but 
of the compression which can be reached by the collid- 
ing system. The next step is to consider now the energy 
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Fig. 13. Excitation function of the ratio R of K + multiplic- 
ities obtained in inclusive Au+Au over C+C reactions. QMD 
[T3U] and IQMD calculations [T3T] are compared to the KaoS 
data 128 . The shaded area indicates thereby the range of un- 
certainty in the theoretical models. In addition IQMD results 
based on an alternative set of elementary K + production cross 
sections are shown. 



dependence of the EOS effect. It is expected to be most 
pronounced most far below threshold because there the 
highest degree of collectivity, reflected in multi-step colli- 
sions, is necessary to overcome the production thresholds. 
The effects become even more evident when the ratio R of 
the kaon multiplicities obtained in Au+Au over C+C re- 
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actions (normalized to the corresponding mass numbers) 
is built |130I128| . Such a ratio has moreover the advan- 
tage that possible uncertainties which might still exist in 
the theoretical calculations should cancel out to large ex- 
tent. This ratio is shown in Fig. f^H Both, soft and hard 
EOS, show an increase of R with decreasing energy. How- 
ever, this increase is much less pronounced when the stiff 
EOS is employed. The strong increase of R can be directly 
related to higher compressible nuclear matter. The com- 
parison to the experimental data from KaoS |128| where 
the increase of R is even more pronounced strongly favors 
a soft equation of state. These findings were confirmed by 
independent IQMD transport calculations of the Nantes 
group [1211. Both, QMD and IQMD included also a re- 
pulsive kaon-nucleon potential as predicted by chiral per- 
turbation theory |129| . The shaded area in the figure can 
be taken as the existing range of uncertainty in the theo- 
retical model description of the considered observable. To 
estimate the stability of the conclusions the IQMD cal- 
culations have been repeated with an alternative set of 
NA: AA i — ► NYK + cross sections 4 which are almost one 
order of magnitude smaller than those used originally, but 
the EOS dependence remained stable |132| . 



6 Constraints from neutron stars 

Measurements of "extreme" values, like large masses or 
radii, huge luminosities etc. as provided by compact stars 
offer good opportunities to gain deeper insight into the 
physics of matter under extreme conditions. There has 
been substantial progress in recent time from the astro- 
physical side. 

The most spectacular observation was probably the 
recent measurements on PSR J0751+1807, a millisecond 
pulsar in a binary system with a helium white dwarf sec- 
ondary, which implies a pulsar mass of 2.1 ± 0.2 (jig's) Mq 
with la (2a) confidence |133| . Therefore a reliable EOS 
has to describe neutron star (NS) masses of at least 1.9 M Q 
(la) in a strong, or 1.6 M (2a) in a weak interpretation. 
This condition limits the softness of EOS in NS matter. 
One might therefore be worried about an apparent con- 
tradiction between the constraints derived from neutron 
stars and those from heavy ion reactions. While heavy 
ion reactions favor a soft EOS, PSR J0751+1807 requires 
a stiff EOS. The corresponding constraints are, however, 
complementary rather than contradictory. Intermediate 
energy heavy ion reactions, e.g. subthreshold kaon pro- 
duction, constrains the EOS at densities up to 2 -=- 3 po 
while the maximum NS mass is more sensitive to the high 
density behaviour of the EOS. Combining the two con- 
straints implies that the EOS should be soft at moderate 
densities and stiff at high densities. Such a behaviour is 
predicted by microscopic many-body calculations (see Fig. 
EJ). DBHF, BHF or variational calculations lead typically 
the maximal NS masses between 2.1 2.3 M Q and are 
therefore in accordance with PSR J0751+1807 |T53| . 

4 Cross section which involve A resonances in the initial of 
final states are not constrained by measurements. 



There exist several other constraints on the nuclear 
EOS which can be derived from observations of compact 
stars, see e.g. |134I135I136| . Among these the one of the 
most promising constraints is the Direct Urea (DU) pro- 
cess which is essentially driven by the proton fraction in- 
side the NS |137| . DU processes, e.g. the neutron /3-decay 
n — > p + e~ + v e , are very efficient regarding their neutrino 
production, even in superfluid NM |138I139| . and cool 
NSs too fast to be in accordance with data from ther- 
mal observable NSs. Therefore one can suppose that no 
DU processes should occur below the upper mass limit 
for "typical" NSs, i.e. M DU > 1.5 M© (1.35 M Q in a 
weak interpretation). These limits come from a popula- 
tion synthesis of young, nearby NSs |140j and masses of 
NS binaries [T35] . 



7 Summary and outlook 

The quest for the nuclear equation of state is one of the 
longstanding problems in physics which has a more than 
50 years history in nuclear structure. Since about 30 years 
one tries to attack this question with heavy ion reactions. 
The exploration of the limits of stability, i.e. the regimes of 
extreme isospin asymmetry, is a relatively new field with 
rapidly growing importance in view of the forthcoming 
generation of radioactive beam facilities. 

The status of theoretical models which make predic- 
tions for the EOS can roughly be summarized as follows: 
phenomcnological density functionals such as Skyrme, Gogny 
or relativistic mean field models provide high precision fits 
to the nuclear chart but extrapolations to supra-normal 
densities or the limits of stability are highly uncertain. 
A more controlled way provide effective field theory ap- 
proaches which became quite popular in recent time. Ef- 
fective chiral field theory allows e.g. a systematic genera- 
tion of two- and many-body nuclear forces. However, these 
approaches are low momentum expansions and when ap- 
plied to the nuclear many-body problem, low density ex- 
pansions. Ab initio calculations for the nuclear many-body 
problem such as variational or Brueckner calculations have 
reached a high degree of sophistication and can serve as 
guidelines for the extrapolation to the regimes of high den- 
sity and/or large isospin asymmetry. Possible future dev- 
elopments are to base such calculations on modern EFT 
potentials and to achieve a more consistent treatment of 
two- and three-body forces. 

If one intends to constrain these models by nuclear re- 
actions one has to account for the reaction dynamics by 
semi-classical transport models of a Boltzmann or molec- 
ular dynamics type. Suitable observables which have been 
found to be sensitive on the nuclear EOS are directed and 
elliptic collective flow pattern and particle production, in 
particular kaon production, at higher energies. Heavy ion 
data suggest that the EOS of symmetric nuclear matter 
shows a soft behavior in the density regime between one to 
about three times nuclear saturation density, which is con- 
sistent with the predictions from many-body calculations. 
Conclusions on the EOS are, however, complicated by the 
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interplay between the density and the momentum depen- 
dence of the nuclear mean field. Data which constrain 
the isospin dependence of the mean field are still scare. 
Promising observables are isospin diffusion, iso-scaling of 
intermediate mass fragments and particle ratios (7r + /7r _ 
and eventually K + /K°). Here the situation will certainly 
improve when the forthcoming radioactive beam facilities 
will be operating. This will also allow to measure the op- 
tical isospin potential in p+A and A+A reactions and 
to obtain more information of the symmetry energy and 
the proton/neutron mass splitting in asymmetric matter. 
From the theoretical side it will be unavoidable to invest 
significant efforts towards the development of quantum 
transport models with consistent off-shell dynamics. 

We would like to thank K. Morawetz, T. Gaitanos and 
M. di Toro for fruitful discussions. 
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